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The free-electron-Uke state observed in a scanning tunneling spectroscopy study of a chiral p(2x4) 
monolayer of glycinate ions on the Cu(lOO) surface [K. Kanazawa et al, J. Am. Chem. Soc. 2007, 
129^ 740] is shown from density functional theory calculations to originate from a Cu Shockley surface 
state at the surface Brillouin zone boundary of the clean surface with highly anisotropic dispersion. 
The presence of the glycinate ions on the surface causes a dramatically enhanced tunneling into this 
surface state that is otherwise not observed in tunneling on the bare surface. 

PACS numbers: 68.37.Ef, 68.43.bc, 73.20.At 



Recently, free-electron-like states arising in layers of 
organic molecules on metal surfaces at energies close to 
the Fermi energy have been observed by scanning tunnel- 
ing and photoemission spectroscopies [U El [3] . Interest 
in such delocalized electronic states is high, because of 
their potential use in molecular electronics applications. 
In systems where the molecule is more or less physisorbed 
the observed free-electron like state was simply identified 
as a metal surface state [3 , whereas in systems where 
the molecule is chemisorbed the origin of the state is not 
properly understood [U [2] . 

In the scanning tunneling spectroscopy (STS) study by 
Kanazawa et al[l\ an anisotropic free-electron-like state 
was observed in a p(2x4) monolayer of glycinate an- 
ions on the Cu(lOO) surface just above the Fermi energy. 
From an analysis of the two-dimensional autocorrelations 
of the differential conductance (^) images , they were 
able to show that the level of dispersion in the [110] and 
[110] directions differs by an order of magnitude. From 
the absence of any observation of a surface state on the 
clean surface in the ^ spectra, they suggested that the 
state originated from an interaction between the glyci- 
nate monolayer and the Cu surface and its anisotropic 
dispersion from the different strength of interaction be- 
tween neighboring molecules in the two directions. There 
is need to carry out electronic structure calculations to 
clarify the origin of this state. 

In this letter we present the results of density func- 
tional theory (DFT) calculations on the p(2x4) glycinate 
monolayer and show that the origin of the anisotropic 
free-electron-like state (AFES) is a Cu Shockley surface 
state (SS). Although Kanazawa et al did not observe a SS 
near the Fermi energy on the bare Cu (100) surface [1] 
a SS has been identified near the Fermi energy at the 
surface Brillouin zone (SBZ) boundary [H [SJ [6l [7] . We 
show that the SS is folded back to the T point in the 
(SBZ) for the p(2x4) structure studied by Kanazawa et 
al^ and that the anisotropic parabolic dispersion of the 
observed AFES is due to enhanced tunneling into the SS, 
mediated by the lowest unoccupied molecular orbitals of 
the glycinate ions. 



Glycine (NII2CH2COOII) is the simplest amino acid 
and its adsorption geometry on the Cu (100) surface is 
well understood following several studies in the last ten 
years with scanning tunneling microscopy (STM) [UIHllQ], 
photoelectron diffraction [10] and DFT calculations pT] . 
The glycine molecule loses a hydrogen atom on adsorp- 
tion becoming a glycinate anion (NIl2CH2COO~) which 
then binds to the Cu surface through the nitrogen atom 
and two oxygen atoms in a tridentate fashion. The 
p(2x4) structure [U |9] has alternating rows of the left 
and right enantiomers of glycinate ions, the rows prop- 
agating along the close-packed rows of Cu atoms (see 
Fig.jl]). In addition to the hetereochiral p(2x4) structure 
a homochiral c(2x4) [8l|9l[T2] structure is also formed at 
the same temperature on the Cu(lOO) surface. 

Density functional theory (DFT) calculations were car- 
ried out using the plane- wave based VASP code [13 . Cal- 
culations were made using the projector augmented wave 
method [14 and the exchange-correlation functional us- 
ing the generalized gradient approximation proposed by 
Perdew and Wang [15 . Glycinate ions were placed on a 
six layer Cu(lOO) slab in a 5.14x10.28x25 A^ supercell. 
A full geometry relaxation was carried out with a plane 
wave cutoff of 400 eV on a 8x4x1 /c-point grid until the 
forces acting on the ions were smaller in magnitude than 
0.01 eV/A. The atoms in the top three Cu layers were 
also allowed to relax. Calculations were also performed 
with a 12 layer bare and adsorbate-covered Cu slab to 
ensure that effects caused by the finite thickness of the 
slab were not important. 

Following the commonly used Tersoff-Hamann approx- 
imation [T6^ the differential conductance measured in 
STS experiments was approximated using the local den- 
sity of states (LDOS) calculated at the position of the 
tip, which was set at 7 A from the top layer of Cu atoms. 
Once the relaxed geometry had been found, the local 
density of states (LDOS) and molecular orbital projected 
density of states (MO-PDOS) was calculated on the orig- 
inal /c-point grid and at 20 points in lines from the f point 
to both the X and Y points in the SBZ (see Fig. [ijc)). 
These calculations used a charge density calculated self- 
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consistently on a 16x8x1 /c-point grid. The MO-PDOS 
was calculated by projecting the density of states of the 
full system onto the wavefunctions for the isolated mono- 
layer in the same supercell at each /c-point. 

Fig. [l] shows the adsorption geometry of the glyci- 
nate molecules following geometry relaxation from one 
close to that found by Mae and Morikawa[lJ] for the 
p(2x4) structure. The resulting Cu-N bond length of 
2.10 A and the Cu-0 bond lengths of 2.08 A and 2.18 A 
are in good agreement with previous theoretical [1T| and 
experimental [10 findings. 
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FIG. 1: Calculated adsorption geometry of glycinate 
molecules on Cu (100) (a) side view and (b) top view. The 
2x4 unit cell is shown from above in white, (c) Labels for 
high symmetry points in the square SBZ of the Cu(lOO) sur- 
face and the rectangular SBZ of the p(2x4) glycinate system. 
The f -X direction corresponds to the direction in real space 
labeled [110] by Kanazawa et a/|l], and F-Y the [110] direc- 
tion. The X^ and Y^ points in the substrate SBZ fold back to 
the F point in the monolayer SBZ. (color online) 



interaction between the glycinate anions and the Cu sur- 
face. Although bands are formed due to interactions 
between neighboring molecules they lie too low in en- 
ergy and do not have the correct dispersion behavior to 
explain free-electron-like state observed by Kanazawa et 
al [1 . 

The calculated LDOS shows a peak which has an 
anisotropic free-electron-like dispersion in very good 
agreement to that found by Kanazawa et al [1] (Figure[3|. 
They found that the dispersion of the AFES giving rise 
to the peak in the auto-correlated ^ spectra was well 
described by parabolic dispersion with effective masses of 
0.061 me and 0.61 me in the [110] and the [110] directions, 
respectively. The dispersion of the peak in the calculated 
LDOS shown in Fig. [3] agrees very well with the experi- 
mentally observed parabolic dispersion, although the ex- 
perimentally observed onset of the peak was at 130 meV 
above the Fermi energy slightly higher than the calcu- 
lated value of 40 meV. The origin of the AFES behind 
the calculated peak can be understood from a scrutiny 
of the electronic states of the bare Cu slab. 




FIG. 2: Bandstructure of a six layer Cu(lOO) slab. The dotted 
curve shows a parabolic dispersion with the effective mass of 
the surface state on the bare Cu surface measured by Kevan[4] 
and the dashed curves with the effective masses observed by 
Kanazawa et al for the state with the glycinate monolayer 
present [T]. 



The electronic structure of molecular monolayer was 
investigated by calculating the MO-PDOS of the ad- 
sorbed system along with the band decomposed DOS 
of an isolated layer of neutral glycinate radicals in the 
same geometry as on the Cu surface. The isolated 
layer is found to have a partially occupied band with a 
width of approximately 0.4 eV across the SBZ formed 
by the overlap of molecular orbitals on neighboring 
molecules. Projection of the DOS of the full adsorbed 
system onto the partially occupied band in the isolated 
glycinate layer shows that the band becomes fully oc- 
cupied in the presence of the Cu surface forming nega- 
tively charged glycinate anions in agreement with previ- 
ous studies [l] [9l [TOl |^. The glycinate states are con- 
siderably broadened upon adsorption reflecting a strong 



Calculations show that the SS of the bare Cu slab has 
a similar energy and anisotropic dispersion to the AFES. 
The SS is clearly observable at the X' point in the cal- 
culated band structure shown in Fig. [2] Two degenerate 
SS are found on each surface of the bare slab; one at the 
X' point and one at the equivalent point for the bare 
surface. The energy splitting in the SS at the X' point 
(see Table |l| due to interactions between the two differ- 
ent sides of the slab is very small even with the six layer 
slab showing that the SS decays quickly into the bulk. 
The calculated dispersion of the SS near the X' point re- 
produces closely the observed parabolic dispersion with 
an effective mass of 0.067±0.01 me in the [110] direction 
measured by Kevan|4j in an angle resolved photoemis- 
sion study. The calculated SS energy of 45 meV below 
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the Fermi energy at the X' point is also in good agree- 
ment with the measured onset of 58±0.05 meV below 
the Fermi energy pi. The effective mass of 0.067±0.01 
me in the [110] direction of the SS is close the measured 
effective mass of 0.061 me of the AFES in the same direc- 
tion. Furthermore, the calculated SS in Fig. [2] shows also 
a parabolic-like dispersion in the [110] direction in close 
agreement with the measured dispersion of the AFES in 
Fig. [2] in the same surface direction with an effective mass 
of 0.61 me- 
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FIG. 3: K-point decomposed LDOS at 7 A above the Cu 
surface in the (a) f-X direction and the (b) f-Y direction. 
The dashed curves are parabolic dispersions with the effective 
masses observed by Kanazawa a/ [T] . The lower set of curves 
are for the fully relaxed adsorbate geometry, the middle set 
with the glycinate ions displaced 0.5 A away from the Cu 
surface and the upper set 1.0 A away. Data from neighboring 
/c-points have been displaced vertically and crosses mark the 
position of the top of the peak at its base. 

The close agreement of the calculated dispersions of 
the SS and the AFES suggests that the observed AFES 
is due to tunneling into the SS. Tunneling into the SS 
must be enhanced by the monolayer, however, since the 
SS cannot be observed on the bare surface using STS [l]. 
Kanazawa et al observed a peak in the ^ spectrum over 
the glycinate monolayer with the steep onset and gradual 
decay with increasing bias characteristic of a dispersive 
state, but no such peak was observed in the ^ spectrum 
over the clean Cu surface. Calculations of the LDOS at 
7 A above the bare Cu slab show no contribution due to 
the SS, only contributions from bulk states around the 
F-point. On the bare substrate the SS at the X' and Y' 



points decay too rapidly into the vacuum region to be 
observed. 

That the calculated AFES arises solely from the molec- 
ular states in the monolayer is ruled out by the behavior 
of the calculated LDOS when changing the distance of 
the glycinate ion monolayer from the Cu surface by 0.5 A 
and 1 A without subsequent relaxation of the geometry. 
The peak with anisotropic dispersion in LDOS (Figure 
[3| is seen to get gradually smaller as the glycinate ions 
are withdrawn from the surface. If the peak were only 
due to tunneling into the Cu surface then the height of 
the peak should stay constant. A peak due to tunnel- 
ing through a purely molecular state would get taller as 
the molecules were brought closer to the tip. This be- 
havior can be seen in the peaks below the Fermi energy. 
The peak with anisotropic dispersion gets smaller as the 
molecules are withdrawn showing that AFES is neither 
a purely molecular state or SS. However, this result sug- 
gests that enhanced tunneling into the SS arises from an 
interaction between the glycinate ions and the Cu sur- 
face. 

An interaction of the SS with the glycinate monolayer 
is confirmed by the fact that the degeneracy of the SS at 
the X' and the Y' points has been lifted and their ener- 
gies are shifted up by 85 meV and 148 meV, respectively, 
upon adsorption. The SS arising from both the X' and 
Y' points on the bare surface are folded back to the f 
point in the 2x4 unit cell. They can be distinguished 
by comparing their dispersion in the f-X (||[110]) and F- 
Y (II [110]) directions. The SS arising from the X' point 
has a smaller effective mass in the F-X direction and that 
arising from the Y' has a smaller effective mass in the F-Y 
direction. The origin of these adsorption-induced energy 
shifts of the SS was investigated by calculating the SS 
energy of the frozen bare slab where Cu atoms were held 
in the geometry they had in the presence the glycinate 
anions. As shown in Table [l[ these energies show that 
although part of the increase in energy of the Shockley 
surface state is due to the adsorbate-induced changes in 
geometry of the Cu surface, part of the energy shift is 
only due to the presence of the glycinate anions. A sub- 
stantial part of the shift in energy must therefore be due 
to direct electronic interaction between the SS and the 
glycinate anions states. The energies of the SS from both 
the X' and Y' points on the bottom side of the slab re- 
mains roughly constant and are degenerate, showing that 
the slab is thick enough to prevent interactions between 
the monolayer and the SS at the bottom of the slab. 

A most interesting result is that only the SS arising 
from the X' point is responsible for the AFES in the 
LDOS in spite of the fact that a larger energy shift is 
found for the SS arising from the Y' point, indicating 
a stronger interaction with the glycinate anions. This 
result explains nicely the observation in the STS exper- 
iments [1 that the AFES has a smaller effective mass 
in the f-X direction than in the f-Y direction. The SS 
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TABLE I: Calculated energies of the Shockley surface state. 





Relaxed 


Frozen 


Glycinate 


Surface, k-point 


bare slab 


bare slab^ 


monolayer 


Top, X' 


-45^ 


+3 


+40 


Top, r 


-45 


-15 


+ 103 


Bottom, X' 


-8 


-5 


-5 


Bottom, 


-8 


-5 


-5 



^Six layer Cu slab with the Cu atoms held in the geometry they 
have when the glycinate atoms are present 
^All energies in meV relative to the Fermi energy of the system. 



arising from the point would have a higher effective 
mass in F-X direction than in the f-Y direction. That 
only the SS from the X' point contribute to the LDOS 
is consistent with the suggestion that the LUMO of the 
glycinate anions plays an important role in the enhanced 
tunneling into the SS. The separate contributions to the 
MO-PDOS from different A:-points in the f -X and f-Y di- 
rections show that the SS giving rise to the AFES in the 
LDOS has non-zero overlap with several of the glycinate 
anion states. The largest overlap is with the first fully 
unoccupied band in the isolated monolayer, the LUMO 
states of the glycinate anions. In contrast, the SS arising 
from the Y' point has no discernible overlap with these 
molecular orbitals. This lack of overlap is not simply a 
symmetry effect. The monolayer has only one non-trivial 
symmetry operation: a combined translation by two Cu 
lattice spacings in the f-Y and a reflection. Because of 
the predominant p character of the SS at the X' and Y' 
points, they are both even under this symmetry opera- 
tion. 
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FIG. 4: Topography of constant orbital density of the AFES 
at the F point with an average distance of 7.0 A from the Cu 
surface. Black (white) represents a height of 6.9 (7.2) A from 
the Cu surface. Each peak is centered over single glycinate 
anion. 



The importance of the glycinate monolayer for the en- 
hanced tunneling into the SS is illustrated by the plot 
in Fig. [4] of the topography of the orbital density of the 



AFES at the F-point in a plane located 7A away from 
the surface. This plot shows that the tunneling into the 
SS is localized to the glycinate ions and centered over the 
protruding hydrogen atom. The localization of the tun- 
neling above the glycinate ions and the overlap of only 
the SS arising from the with X' point with the LUMO 
suggests that this small overlap is sufficient for resonant 
tunneling through the LUMO into the SS beneath the 
molecules. 

In conclusion, we have shown from our density func- 
tional calculations that an organic molecular monolayer 
can dramatically enhance the tunneling into a metal sur- 
face state that is otherwise not observed in scanning tun- 
neling spectroscopy (STS). This finding is of importance 
in understanding the origin of free-electron like states in 
tunneling through molecular monolayers. 
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